LA-GEM: imputation of gene expression with incorporation of Local Ancestry.

Gene imputation and TWAS have become a staple in the genomics medicine discovery space; helping to identify genes whose regulation effects may contribute to disease susceptibility. However, the cohorts on which these methods are built are overwhelmingly of European Ancestry. This means that the unique regulatory variation that exist in non-European populations, specifically African Ancestry populations, may not be included in the current models. Moreover, African Americans are an admixed population, with a mix of European and African segments within their genome. No gene imputation model thus far has incorporated the effect of local ancestry (LA) on gene expression imputation. As such, we created LA-GEM which was trained and tested on a cohort of 60 African American hepatocyte primary cultures. Uniquely, LA-GEM include local ancestry inference in its prediction of gene expression. We compared the performance of LA-GEM to PrediXcan trained the same dataset (with no inclusion of local ancestry) We were able to reliably predict the expression of 2559 genes (1326 in LA-GEM and 1236 in PrediXcan). Of these, 546 genes were unique to LA-GEM, including the CYP3A5 gene which is critical to drug metabolism. We conducted TWAS analysis on two African American clinical cohorts with pharmacogenomics phenotypic information to identity novel gene associations. In our IWPC warfarin cohort, we identified 17 transcriptome-wide significant hits. No gene reached are prespecified significance level in the clopidogrel cohort. We did see suggestive association with RAS3A to P2RY12 Reactivity Units (PRU), a clinical measure of response to anti-platelet therapy. This method demonstrated the need for the incorporation of LA into study in admixed populations.

Investigation of genomic and transcriptomic risk factors in clopidogrel response in African Americans.

Clopidogrel, an anti-platelet drug, used to prevent thrombosis after percutaneous coronary intervention. Clopidogrel resistance results in recurring ischemic episodes, with African Americans suffering disproportionately. The aim of this study was to identify biomarkers of clopidogrel resistance in African American patients. We conducted a genome-wide association study, including local ancestry adjustment, in 141 African Americans on clopidogrel to identify associations with high on-treatment platelet reactivity (HTPR). We validated genome-wide and suggestive hits in an independent cohort of African American clopidogrel patients (N = 823) from the Million Veteran's Program (MVP) along with in vitro functional follow up. We performed differential gene expression (DGE) analysis in whole blood with functional follow-up in MEG-01 cells. We identified rs7807369, within thrombospondin 7A (THSD7A), as significantly associated with increasing risk of HTPR (p = 4.56 × 10-9). Higher THSD7A expression was associated with HTPR in an independent gene expression cohort of clopidogrel treated patients (p = 0.004) and supported by increased gene expression on THSD7A in primary human endothelial cells carrying the risk haplotype. Two SNPs (rs1149515 and rs191786) were validated in the MVP cohort. DGE analysis identified an association with decreased LAIR1 expression to HTPR. LAIR1 knockdown in a MEG-01 cells resulted in increased expression of SYK and AKT1, suggesting an inhibitory role of LAIR1 in the Glycoprotein VI pathway. Notably, the CYP2C19 variants showed no association with clopidogrel response in the discovery or MVP cohorts. In summary, these finding suggest that other variants outside of CYP2C19 star alleles play an important role in clopidogrel response in African Americans.

Discovery of Ancestry-specific Variants Associated with Clopidogrel Response among Caribbean Hispanics.

Background: High on-treatment platelet reactivity (HTPR) with clopidogrel is predictive of ischemic events in adults with coronary artery disease. Despite strong data suggesting HTPR varies with ethnicity, including clinical and genetic variables, no genome-wide association study (GWAS) of clopidogrel response has been performed among Caribbean Hispanics. This study aimed to identify genetic predictors of HTPR in a cohort of Caribbean Hispanic cardiovascular patients from Puerto Rico. Methods: Local Ancestry inference (LAI) and traditional GWASs were performed on a cohort of 511 clopidogrel-treated patients, stratified based on their P2Y12 reaction units (PRU) into responders and non-responders (HTPR). Results: The LAI GWAS identified variants within the CYP2C19 region associated with HTPR, predominantly driven by individuals of European ancestry and absent in those with native ancestry. Incorporating local ancestry adjustment notably enhanced our ability to detect associations. While no loci reached traditional GWAS significance, three variants showed suggestive significance at chromosomes 3, 14 and 22 (OSBPL10 rs1376606, DERL3 rs5030613, and RGS6 rs9323567). In addition, a variant in the UNC5C gene on chromosome 4 was associated with an increased risk of HTPR. These findings were not identified in other cohorts, highlighting the unique genetic landscape of Caribbean Hispanics. Conclusion: This is the first GWAS of clopidogrel response in Hispanics, confirming the relevance of the CYP2C19 cluster, particularly among those with European ancestry, and also identifying novel markers in a diverse patient population. Further studies are warranted to replicate our findings in other diverse cohorts and meta-analyses.

Molecular Identification and Evolutionary Divergence of the Sri Lankan Sambar Deer, Rusa unicolor (Kerr 1792).

The Sambar is one of the largest deer species distributed mainly in Asia, and it has been listed as a vulnerable species. Taxonomy based on morphological characterization has been the gold standard method used to identify the Sambar deer species. Yet, morphological identification is challenging and requires expertise. To conduct species identification and taxonomic decisions, we performed the molecular identification of R. unicolor found in Sri Lanka using DNA barcodes, COI, and Cyt b to compare the Sri Lankan R. unicolor with the Indian R. unicolor and other R. unicolor subspecies. We obtained mitochondrial DNA sequences from COI and Cyt b from blood samples collected from the wet zone in Sri Lanka. A phylogenetic tree was constructed based on the Bayesian analyses using MrBayes 3.2.7. Molecular dating was implemented in Bayesian Evolutionary Analysis Sampling Trees (BEAST v1.8.2) on the concatenated sequence using a log-normal relaxed clock and Yule species tree prior, with four categories. The results showed that the Sri Lankan R. unicolor is genetically different from the Indian R. unicolor and other R. unicolor subspecies. The divergence occurred approximately 1.1 MYA (million years ago) in the Pleistocene era. The results are essential for designing new conservation platforms for these Sambar deer species.

Prostate cancer cell-platelet bidirectional signaling promotes calcium mobilization, invasion and apoptotic resistance via distinct receptor-ligand pairs.

Platelets play a crucial role in cancer and thrombosis. However, the receptor-ligand repertoire mediating prostate cancer (PCa) cell-platelet interactions and ensuing consequences have not been fully elucidated. Microvilli emanating from the plasma membrane of PCa cell lines (RC77 T/E, MDA PCa 2b) directly contacted individual platelets and platelet aggregates. PCa cell-platelet interactions were associated with calcium mobilization in platelets, and translocation of P-selectin and integrin αIIbß3 onto the platelet surface. PCa cell-platelet interactions reciprocally promoted PCa cell invasion and apoptotic resistance, and these events were insensitive to androgen receptor blockade by bicalutamide. PCa cells were exceedingly sensitive to activation by platelets in vitro, occurring at a PCa cell:platelet coculture ratio as low as 1:10 (whereas PCa patient blood contains 1:2,000,000 per ml). Conditioned medium from cocultures stimulated PCa cell invasion but not apoptotic resistance nor platelet aggregation. Candidate transmembrane signaling proteins responsible for PCa cell-platelet oncogenic events were identified by RNA-Seq and broadly divided into 4 major categories: (1) integrin-ligand, (2) EPH receptor-ephrin, (3) immune checkpoint receptor-ligand, and (4) miscellaneous receptor-ligand interactions. Based on antibody neutralization and small molecule inhibitor assays, PCa cell-stimulated calcium mobilization in platelets was found to be mediated by a fibronectin1 (FN1)-αIIbß3 signaling axis. Platelet-stimulated PCa cell invasion was facilitated by a CD55-adhesion G protein coupled receptor E5 (ADGRE5) axis, with contribution from platelet cytokines CCL3L1 and IL32. Platelet-stimulated PCa cell apoptotic resistance relied on ephrin-EPH receptor and lysophosphatidic acid (LPA)-LPA receptor (LPAR) signaling. Of participating signaling partners, FN1 and LPAR3 overexpression was observed in PCa specimens compared to normal prostate, while high expression of CCR1 (CCL3L1 receptor), EPHA1 and LPAR5 in PCa was associated with poor patient survival. These findings emphasize that non-overlapping receptor-ligand pairs participate in oncogenesis and thrombosis, highlighting the complexity of any contemplated clinical intervention strategy.

Implementation of pharmacogenomics into inpatient general medicine.

Pharmacogenomics is a crucial piece of personalized medicine. Preemptive pharmacogenomic testing is only used sparsely in the inpatient setting and there are few models to date for fostering the adoption of pharmacogenomic treatment in the inpatient setting. We created a multi-institutional project in Chicago to enable the translation of pharmacogenomics into inpatient practice. We are reporting our implementation process and barriers we encountered with solutions. This study, 'Implementation of Point-of-Care Pharmacogenomic Decision Support Accounting for Minority Disparities', sought to implement pharmacogenomics into inpatient practice at three sites: The University of Chicago, Northwestern Memorial Hospital, and the University of Illinois at Chicago. This study involved enrolling African American adult patients for preemptive genotyping across a panel of actionable germline variants predicting drug response or toxicity risk. We report our approach to implementation and the barriers we encountered engaging hospitalists and general medical providers in the inpatient pharmacogenomic intervention. Our strategies included: a streamlined delivery system for pharmacogenomic information, attendance at hospital medicine section meetings, use of physician and pharmacist champions, focus on hospitalists' care and optimizing system function to fit their workflow, hand-offs, and dealing with hospitalists turnover. Our work provides insights into strategies for the initial engagement of inpatient general medicine providers that we hope will benefit other institutions seeking to implement pharmacogenomics in the inpatient setting.

Leveraging drug perturbation to reveal genetic regulators of hepatic gene expression in African Americans.

Expression quantitative locus (eQTL) studies have paved the way in identifying genetic variation impacting gene expression levels. African Americans (AAs) are disproportionately underrepresented in eQTL studies, resulting in a lack of power to identify population-specific regulatory variants especially related to drug response. Specific drugs are known to affect the biosynthesis of drug metabolism enzymes as well as other genes. We used drug perturbation in cultured primary hepatocytes derived from AAs to determine the effect of drug treatment on eQTL mapping and to identify the drug response eQTLs (reQTLs) that show altered effect size following drug treatment. Whole-genome genotyping (Illumina MEGA array) and RNA sequencing were performed on 60 primary hepatocyte cultures after treatment with six drugs (Rifampin, Phenytoin, Carbamazepine, Dexamethasone, Phenobarbital, and Omeprazole) and at baseline (no treatment). eQTLs were mapped by treatment and jointly with Meta-Tissue. We found varying transcriptional changes across different drug treatments and identified Nrf2 as a potential general transcriptional regulator. We jointly mapped eQTLs with gene expression data across all drug treatments and baseline, which increased our power to detect eQTLs by 2.7-fold. We also identified 2,988 reQTLs (eQTLs with altered effect size after drug treatment). reQTLs were more likely to overlap transcription factor binding sites, and we uncovered reQTLs for drug metabolizing genes such as CYP3A5. Our results provide insights into the genetic regulation of gene expression in hepatocytes through drug perturbation and provide insight into SNPs that effect the liver's ability to respond to transcription upregulation.

Genomewide Association Study Identifies Copy Number Variants Associated With Warfarin Dose Response and Risk of Venous Thromboembolism in African Americans.

The anticoagulant warfarin is commonly used to control and prevent thrombotic disorders, such as venous thromboembolism (VTE), which disproportionately afflicts African Americans. Despite the importance of copy number variants (CNVs), few studies have focused on characterizing and understanding their role in drug response and disease risk among African Americans. In this study, we conduct the first genome-wide analysis of CNVs to more comprehensively account for the contribution of genetic variation in warfarin dose requirement and VTE risk among African Americans. We used hidden Markov models to detect CNVs from high-throughput single-nucleotide polymorphism arrays for 340 African American participants in the International Warfarin Pharmacogenetics Consortium. We identified 11,570 CNVs resulting in 2,038 copy number variable regions (CNVRs) and found 3 CNVRs associated with warfarin dose requirement and 3 CNVRs associated with VTE risk in African Americans. CNVRs 1q31.2del and 6q14.1del were associated with increased warfarin dose requirement (ß = 11.18 and 4.94, respectively; Pemp  = < 0.002); CNVR 19p13.31del was associated with decreased warfarin dose requirement (ß = -1.41, Pemp  = 0.0004); CNVRs (2p22.1del and 5q35.1-q35.2del) were found to be associated with increased risk of VTE (odds ratios (ORs) = 1.88 and 14.9, respectively; Pemp ≤0.02); and CNVR 10q26.12del was associated with a decreased risk of VTE (OR = 0.6; Pemp  = 0.05). Modeling of the 10q26.12del in HepG2 cells revealed that this deletion results in decreased fibrinogen gene expression, decreased fibrinogen and WDR11 protein levels, and decreased secretion of fibrinogen into the extracellular matrix. We found robust evidence that CNVRs could contribute to warfarin dose requirement and risk of VTE in African Americans and for 10q26.3del describe a possible pathogenic mechanism.

The Role of Global and Local Ancestry on Clopidogrel Response in African Americans.

Pharmacogenomics has long lacked dedicated studies in African Americans, resulting in a lack of indepth data in this populations. The ACCOuNT consortium has collected a cohort of 167 African American patients on steady state clopidogrel with the goal of discovering population specific variation that may contribute to the response of this anti-platelet agent. Here we analyze the role of both global and local ancestry on the clinical phenotypes of P2Y12 reaction units (PRU) and high on-treatment platelet reactivity (HTPR) in this cohort. We found that local ancestry at the TSS of three genes, IRS-1, ABCB1 and KDR were nominally associated with PRU, and local ancestry-adjusted SNP association identified variants in ITGA2 associated to increased PRU. These finding help to explain the variability in drug response seen in African Americans, especially as few studies on genes outside of CYP2C19 has been conducted in this population.

Local Ancestry-Informed Candidate Pathway Analysis of Warfarin Stable Dose in Latino Populations.

Accuracy of warfarin dose prediction algorithms may be improved by including data from diverse populations in genetic studies of dose variability. Here, we surveyed single nucleotide polymorphisms in vitamin K-related genetic pathways for association with warfarin dose requirements in two admixed Latino populations in standard-principal component adjusted and contemporary-local ancestry adjusted regression models. A total of five variants from vitamin K-related genes/pathways were associated with warfarin dose in both cohorts (P < 0.0125) in standard models. Local ancestry-adjusted analysis unveiled 35 associated variants with absolute effects ranging from ß = 9.04 ( ±2.23) to 39.18 ( ±10.89) per ancestral allele in the discovery cohort and ß = 6.47 (± 2.02) to 17.82 (± 6.83) in the replication cohort. Importantly, we demonstrate the technical validity of the Tractor model in cohorts with admixed ancestry from three founder populations and bring attention to the technical hurdles obstructing the inclusion of diverse, especially admixed, populations in pharmacogenomic research.

Multi-Omics Studies in Historically Excluded Populations: The Road to Equity.

Over the past few decades, genomewide association studies (GWASs) have identified the specific genetics variants contributing to many complex diseases by testing millions of genetic variations across the human genome against a variety of phenotypes. However, GWASs are limited in their ability to uncover mechanistic insight given that most significant associations are found in non-coding region of the genome. Furthermore, the lack of diversity in studies has stymied the advance of precision medicine for many historically excluded populations. In this review, we summarize most popular multi-omics approaches (genomics, transcriptomics, proteomics, and metabolomics) related to precision medicine and highlight if diverse populations have been included and how their findings have advance biological understanding of disease and drug response. New methods that incorporate local ancestry have been to improve the power of GWASs for admixed populations (such as African Americans and Latinx). Because most signals from GWAS are in the non-coding region, other machine learning and omics approaches have been developed to identify the potential causative single-nucleotide polymorphisms and genes that explain these phenotypes. These include polygenic risk scores, expression quantitative trait locus mapping, and transcriptome-wide association studies. Analogous protein methods, such as proteins quantitative trait locus mapping, proteome-wide association studies, and metabolomic approaches provide insight into the consequences of genetic variation on protein abundance. Whereas, integrated multi-omics studies have improved our understanding of the mechanisms for genetic association, we still lack the datasets and cohorts for historically excluded populations to provide equity in precision medicine and pharmacogenomics.

Warfarin Dosing in Patients With CYP2C9*5 Variant Alleles.

Pharmacogenetic dosing improves the accuracy of warfarin dosing, but current pharmacogenetic dosing algorithms are less accurate in populations of African ancestry. The cytochrome P450 2C9*5 (CYP2C9*5) allele is found almost exclusively in populations of African ancestry, and in vitro studies suggest CYP2C9*5 is associated with reduced clearance of warfarin. The clinical relevance of this single-nucleotide variation (SNV) (formerly SNP) is uncertain. In this multicentered study of 2,298 patients (49% female, 35% Black) taking warfarin, we quantified the association between the CYP2C9*5 allele and warfarin requirements. The CYP2C9*5 SNV was present in 2.3% of Black and 0.07% of White patients. Without taking CYP2C9*5 into account, pharmacogenetic algorithms that include other SNVs overestimated the warfarin dose by 30% (95% confidence interval (19-40%), P < 0.001), an average of 1.87 mg/day (SD 1.64) in heterozygotes (P < 0.001). Noncarriers required a slightly (0.23 mg/day, SD 2.09) higher than predicted dose. Genotyping for CYP2C9*5 corrected the potential overdose and halved overall dosing error in heterozygotes. Patients carrying CYP2C9*5 require a clinically relevant reduction in warfarin dose. Given the potential to improve the accuracy and safety of warfarin dosing in populations of African ancestry, we have incorporated this SNV into a nonprofit website to assist warfarin initiation (www.WarfarinDosing.org).

Apixaban Concentrations in Routine Clinical Care of Older Adults With Nonvalvular Atrial Fibrillation.

BACKGROUND: Direct-acting oral anticoagulants are first-line agents for prevention of stroke in patients with nonvalvular atrial fibrillation (NVAF), but data are limited for the oldest patients, and with reduced dosing. OBJECTIVES: To determine steady-state apixaban peak and trough concentrations during routine care of older adults with NVAF, compare concentrations to clinical trial concentrations, and explore factors associated with concentrations. METHODS: A cross-sectional study of medically stable older adults with NVAF (≥75 years or ≥70 years if Black) receiving apixaban. Peak (2-4.4 hours post-dose) and trough (before next dose) concentrations were determined by anti-Xa activity calibrated chromogenic assay. Patient characteristics associated with concentrations were determined by multivariate modeling. RESULTS: The median age of patients (n = 115) was 80 (interquartile range: 77-84) years. The cohort comprised 46 women and 69 men; of which 98 are White, 11 Black, and 6 Asian. With 5 mg twice daily per labelling (n = 88), peak concentrations were higher in women: 248 ± 105 vs 174 ± 67 ng/mL in men (P < 0.001) and exceeded expected 95% range in 6 of 30 vs 0 of 55 men (P = 0.002). With 2.5 mg twice daily per label (n = 11), concentrations were <5 mg twice daily (peak: 136 ± 87 vs 201 ± 90 ng/mL, P = 0.026; trough: 65 ± 28 vs 109 ± 56 ng/mL, P < 0.001), but not different than 2.5 mg twice daily without reduction criteria (n = 13; peak: 132 ± 88; trough: 65 ± 31 ng/mL). Covariates associated with concentrations included sex, number of daily medications, and creatinine clearance. CONCLUSIONS: Older women had higher than expected peak apixaban concentrations, and 2.5 mg twice daily produced lower concentrations than standard dosing. Factors not currently included in dosing recommendations affected concentrations. The impact of apixaban concentrations on outcomes needs evaluation.

Genetic association of primary nonresponse to anti-TNFα therapy in patients with inflammatory bowel disease.

OBJECTIVES: Primary nonresponse (PNR) to antitumor necrosis factor-α (TNFα) biologics is a serious concern in patients with inflammatory bowel disease (IBD). We aimed to identify the genetic variants associated with PNR. PATIENTS AND METHODS: Patients were recruited from outpatient GI clinics and PNR was determined using both clinical and endoscopic findings. A case-control genome-wide association study was performed in 589 IBD patients and associations were replicated in an independent cohort of 293 patients. Effect of the associated variant on gene expression and TNFα secretion was assessed by cell-based assays. Pleiotropic effects were investigated by Phenome-wide association study (PheWAS). RESULTS: We identified rs34767465 as associated with PNR to anti-TNFα therapy (odds ratio: 2.07, 95% CI, 1.46-2.94, P = 2.43 × 10-7, [replication odds ratio: 1.8, 95% CI, 1.04-3.16, P = 0.03]). rs34767465 is a multiple-tissue expression quantitative trait loci for FAM114A2. Using RNA-sequencing and protein quantification from HapMap lymphoblastoid cell lines (LCLs), we found a significant decrease in FAM114A2 mRNA and protein expression in both heterozygous and homozygous genotypes when compared to wild type LCLs. TNFα secretion was significantly higher in THP-1 cells [differentiated into macrophages] with FAM114A2 knockdown versus controls. Immunoblotting experiments showed that depletion of FAM114A2 impaired autophagy-related pathway genes suggesting autophagy-mediated TNFα secretion as a potential mechanism. PheWAS showed rs34767465 was associated with comorbid conditions found in IBD patients (derangement of joints [P = 3.7 × 10-4], pigmentary iris degeneration [P = 5.9 × 10-4], diverticulum of esophagus [P = 7 × 10-4]). CONCLUSIONS: We identified a variant rs34767465 associated with PNR to anti-TNFα biologics, which increases TNFα secretion through mechanism related to autophagy. rs34767465 may also explain the comorbidities associated with IBD.

Applicability of Pharmacogenomically Guided Medication Treatment during Hospitalization of At-Risk Minority Patients.

Known disparities exist in the availability of pharmacogenomic information for minority populations, amplifying uncertainty around clinical utility for these groups. We conducted a multi-site inpatient pharmacogenomic implementation program among self-identified African-Americans (AA; n = 135) with numerous rehospitalizations (n = 341) from 2017 to 2020 (NIH-funded ACCOuNT project/clinicaltrials.gov#NCT03225820). We evaluated the point-of-care availability of patient pharmacogenomic results to healthcare providers via an electronic clinical decision support tool. Among newly added medications during hospitalizations and at discharge, we examined the most frequently utilized medications with associated pharmacogenomic results. The population was predominantly female (61%) with a mean age of 53 years (range 19-86). On average, six medications were newly prescribed during each individual hospital admission. For 48% of all hospitalizations, clinical pharmacogenomic information was applicable to at least one newly prescribed medication. Most results indicated genomic favorability, although nearly 29% of newly prescribed medications indicated increased genomic caution (increase in toxicity risk/suboptimal response). More than one of every five medications prescribed to AA patients at hospital discharge were associated with cautionary pharmacogenomic results (most commonly pantoprazole/suboptimal antacid effect). Notably, high-risk pharmacogenomic results (genomic contraindication) were exceedingly rare. We conclude that the applicability of pharmacogenomic information during hospitalizations for vulnerable populations at-risk for experiencing health disparities is substantial and warrants continued prospective investigation.

Differences in the Platelet mRNA Landscape Portend Racial Disparities in Platelet Function and Suggest Novel Therapeutic Targets.

The African American (AA) population displays a 1.6 to 3-fold higher incidence of thrombosis and stroke mortality compared with European Americans (EAs). Current antiplatelet therapies target the ADP-mediated signaling pathway, which displays significant pharmacogenetic variation for platelet reactivity. The focus of this study was to define underlying population differences in platelet function in an effort to identify novel molecular targets for future antiplatelet therapy. We performed deep coverage RNA-Seq to compare gene expression levels in platelets derived from a cohort of healthy volunteers defined by ancestry determination. We identified > 13,000 expressed platelet genes of which 480 were significantly differentially expressed genes (DEGs) between AAs and EAs. DEGs encoding proteins known or predicted to modulate platelet aggregation, morphology, or platelet count were upregulated in AA platelets. Numerous G-protein coupled receptors, ion channels, and pro-inflammatory cytokines not previously associated with platelet function were likewise differentially expressed. Many of the signaling proteins represent potential pharmacologic targets of intervention. Notably, we confirmed the differential expression of cytokines IL32 and PROK2 in an independent cohort by quantitative real-time polymerase chain reaction, and provide functional validation of the opposing actions of these two cytokines on collagen-induced AA platelet aggregation. Using Genotype-Tissue Expression whole blood data, we identified 516 expression quantitative trait locuses with Fst values > 0.25, suggesting that population-differentiated alleles may contribute to differences in gene expression. This study identifies gene expression differences at the population level that may affect platelet function and serve as potential biomarkers to identify cardiovascular disease risk. Additionally, our analysis uncovers candidate novel druggable targets for future antiplatelet therapies.

Incorporation of DNA methylation into eQTL mapping in African Americans.

Epigenetics is a reversible molecular mechanism that plays a critical role in many developmental, adaptive, and disease processes. DNA methylation has been shown to regulate gene expression and the advent of high throughput technologies has made genome-wide DNA methylation analysis possible. We investigated the effect of DNA methylation on eQTL mapping (methylation-adjusted eQTLs), by incorporating DNA methylation as a SNP-based covariate in eQTL mapping in African American derived hepatocytes. We found that the addition of DNA methylation uncovered new eQTLs and eGenes. Previously discovered eQTLs were significantly altered by the addition of DNA methylation data suggesting that methylation may modulate the association of SNPs to gene expression. We found that methylation-adjusted eQTLs that were less significant compared to PC-adjusted eQTLs were enriched in lipoprotein measurements (FDR=0.0040), immune system disorders (FDR = 0.0042), and liver enzyme measurements (FDR=0.047), suggesting that DNA methylation modulates the genetic regulation of these phenotypes. Our methylation-adjusted eQTL analysis also uncovered novel SNP-gene pairs. For example, we found that the SNP, rs1332018, was associated to GSTM3. GSTM3 expression has been linked to Hepatitis B which African Americans suffer from disproportionately. Our methylation-adjusted method adds new understanding to the genetic basis of complex diseases that disproportionally affect African Americans.

Discovery of novel hepatocyte eQTLs in African Americans.

African Americans (AAs) are disproportionately affected by metabolic diseases and adverse drug events, with limited publicly available genomic and transcriptomic data to advance the knowledge of the molecular underpinnings or genetic associations to these diseases or drug response phenotypes. To fill this gap, we obtained 60 primary hepatocyte cultures from AA liver donors for genome-wide mapping of expression quantitative trait loci (eQTL) using LAMatrix. We identified 277 eGenes and 19,770 eQTLs, of which 67 eGenes and 7,415 eQTLs are not observed in the Genotype-Tissue Expression Project (GTEx) liver eQTL analysis. Of the eGenes found in GTEx only 25 share the same lead eQTL. These AA-specific eQTLs are less correlated to GTEx eQTLs. in effect sizes and have larger Fst values compared to eQTLs found in both cohorts (overlapping eQTLs). We assessed the overlap between GWAS variants and their tagging variants with AA hepatocyte eQTLs and demonstrated that AA hepatocyte eQTLs can decrease the number of potential causal variants at GWAS loci. Additionally, we identified 75,002 exon QTLs of which 48.8% are not eQTLs in AA hepatocytes. Our analysis provides the first comprehensive characterization of AA hepatocyte eQTLs and highlights the unique discoveries that are made possible due to the increased genetic diversity within the African ancestry genome.

Why African Americans say "No": A Study of Pharmacogenomic Research Participation.

Objective: To identify reasons for nonparticipation by African Americans in cardiovascular pharmacogenomic research. Design: Prospective, open-ended, qualitative survey. Setting: Research staff approached patients eligible for the Discovery Project of The African American Cardiovascular pharmacogenomics CONsorTium in the inpatient or outpatient setting at four different institutions during September and October 2018. Participants: Potential Discovery Project participants self-identified as African American, aged >18 years, were on one of five cardiovascular drugs of interest, and declined enrollment in the Discovery Project. Main Outcome Measures: Reasons for nonparticipation. Methods: After declining participation in the Discovery Project, patients were asked, "What are your reasons for not participating?" We analyzed their responses using a directed content analytic approach. Ultimately, responses were coded into one of nine categories and analyzed using descriptive statistics. Results: Of the 194 people approached for the Discovery Project during an eight-week period, 82 declined participation and provided information for this study. The most common reason for refusal was concern about the amount of blood drawn (19.5%). The next most common reasons for refusal to participate included concerns about genetic testing (14.6%) and mistrust of research (12.2%). Across study sites, significantly more patients enrolled in the inpatient than outpatient setting (P<.001). Significantly more women and younger individuals declined participation due to concerns about genetic testing and too little compensation (P<.05). Conclusions: Collection of blood samples and concerns about genetic testing are obstacles for the recruitment of African Americans to pharmacogenomics studies. Efforts to overcome these barriers to participation are needed to improve representation of minorities in pharmacogenomic research. Enrolling participants from inpatient populations may be a solution to bolster recruitment efforts.

Impact of CYP2C9-Interacting Drugs on Warfarin Pharmacogenomics.

Precise dosing of warfarin is important to achieve therapeutic benefit without adverse effects. Pharmacogenomics explains some interindividual variability in warfarin response, but less attention has been paid to drug-drug interactions in the context of genetic factors. We investigated retrospectively the combined effects of cytochrome P450 (CYP)2C9 and vitamin K epoxide reductase complex (VKORC)1 genotypes and concurrent exposure to CYP2C9-interacting drugs on long-term measures of warfarin anticoagulation. Study participants predicted to be sensitive responders to warfarin based on CYP2C9 and VKORC1 genotypes, had significantly greater international normalized ratio (INR) variability over time. Participants who were concurrently taking CYP2C9-interacting drugs were found to have greater INR variability and lesser time in therapeutic range. The associations of INR variability with genotype were driven by the subgroup not exposed to interacting drugs, whereas the effect of interacting drug exposure was driven by the subgroup categorized as normal responders. Our findings emphasize the importance of considering drug interactions in pharmacogenomic studies.